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Open quantum systems have attracted great attention, since inevitable coupling between 
quantum systems and their environment greatly affects the features of interest of these 
systems. Quantum discord, is a measure of the total nonclassical correlation in a quan- 
tum system that includes, but is not exclusive to, the distinct property of quantum 
entanglement. Quantum discord can exist in separated quantum states and plays an im- 
portant role in many fundamental physics problems and practical quantum information 
tasks. There have been numerous investigations on quantum discord and its counterpart 
classical correlation. This short review focuses on highlighting the system-environment 
dynamics of two-qubit quantum discord and the influence of initial system-environment 
correlations on the dynamics of open quantum systems. The external control effect on 
the dynamics of open quantum systems are involved. Several related experimental works 
are discussed. 

Keywords: quantum discord; classical correlation; system-environment coupling. 



1. Introduction 

Interest in open quantum systems arises naturally from the fact that any quantum 
system is always unavoidably coupled to its surrounding environment^. This makes 
investigation of the preparation, processing, and application of features in open 
quantum systems both fundamentally and practically important. One of the most 
remarkable properties of quantum systems is the existence of correlations that do not 
have corresponding classical counterparts. Entanglement, a special kind of nonclas- 
sical correlation, is one such that has been widely studiec03 and found to be a useful 
resource in quantum communications and quantum computation^. Entanglement is 
fragile, and significant effort has gone into investigating its behavior in open quan- 
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turn systems^. One unusual dynamic entanglement behavior, entanglement sudden 
death, occurs when entanglement completely disappears at a finite point in the sys- 
tem's evolutiorP. Recent studies have discovered quantum correlations other than 
entanglement, with some of these other nonclassical correlations even existing in 
separable quantum states. Quantum correlations without entanglement could po- 
tentially play an important role in implementing quantum information tasks that 
benefit from quantum advantages such as information locking-^, DQC1 (determin- 
istic quantum computation with one quantum bit) algorithm applications^!, quan- 
tum state discriminatiorllHliU, quantum metrologjff^, and remote quantum state 
operations^. Quantum discord^ has also attracted significant attention as a way 
to quantify these total quantum correlations encoded in a quantum systeirli£li£l 

This short review focuses on the dynamics of two-qubit quantum discord un- 
der system-environment coupling, dividing the work into five sections. We first 
introduce the original definitions of quantum discord and classical correlation in 
Part 2, as well as introducing other popular related nonclassical measures. In Part 
3, we review recent theoretical developments of two-qubit quantum discord under 
system-environment coupling, discussing the system's dynamic properties and ex- 
ternal control effects. In Part 4, several related experiments are reviewed in detail, 
followed by a conclusion in Part 5. 

2. Quantum discord and classical correlation 

In classical information theory, the correlations between two different random vari- 
ables A and B, with probability distributions p t and corresponding to the obtained 
values of Oj and 6j, can be characterized by mutual informational 

I(A:B)=H(A)+H(B)-H(A,B), (1) 

In Eq. ([IJ H(A) (H(B)) is the Shannon entropy representing the uncertainty of 
the outcomes of A (B), such that H(A) = — E i p i log 2 Pi (H(B) = —T, i q i \og 2 q i ). 
Note that H(A,B) = —"Sijpij log 2 Pij is the joint entropy where pij represents the 
joint probability of the outcomes of A and B being ai and bj, respectively, and 
Pi = 'SjPij {qj — TiiPij). According to the Bayes rule, the conditional probability 
can be written as p^j — pij/qj. As a result, classical mutual information has the 
equivalent form 

J(A:B)=H(A)-H{A\B), (2) 

where H(A\B) is the conditional entropy of A when B is known, and H(A\B) = 
-EijPij log 2Pilj = H{A, B) ~ H(B). 

In the quantum case, the density matrix p fully characterizes the state of a 
system. The information of a certain state is given by the von Neumann entropy 
S(p), where S{p) = — Tr(/olog 2 p). For a bipartite quantum system with subsystems 
A and B, the quantum mutual information is represented as 



Z(pab) = S(p A ) + S{p B ) - S(p A b), 



(3) 
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where pa = ^bPab (pb = TtaPab) is the reduced density matrix of partition A 
(B). Equation ^ is a direct quantum extension of Eq. (JlJ. Because a measure- 
ment generally disturbs the quantum state, the conditional entropy of subsystem A 
changes when we choose different measurements on subsystem B, and the extension 
of Eq. ([2| to the quantum region is not so direct. 

Consider a set of measurements {M B } on subsystem B. The reduced state of 
subsystem A conditioned on the measurement labeled by k becomes 

p\ = —Tr B [(l A ® M%) PAB {1 A ® M k B )\ (4) 

with probability r k = TiabYk^-a ® M b )pab(1a <8 M B ) t] and with 1 A representing 
the identical operator on the subsystem A. As a result, the conditional entropy of 
subsystem A due to the measurement on B can be defined as 

S(A\B) = Z k r k S(p k A ). (5) 

The quantum generalization of Eq. ([2| becomes 

J(pab) =S( Pa )-S(A\B). (6) 

and the maximum extractable information from the measurement on B 

C(pab) = m&x[J( P AB)} = m&x[S{p A ) - S(A\B)}. (7) 

M B M B 

is then defined as the classical correlation of pab^- The two generalized quantum 
mutual information sets are generally not equal to each other and their difference 
is defined as the famous quantum discord^ 

Q{ P ab) =1{pab)-C{ P ab)- (8) 

The quantum mutual information I(pab) is equal to the sum of classical correlation 
C{pab) and quantum discord Q(pab), an d which is typically used to quantify the 
total correlations for bipartite quantum systems r iH122l 

From Eq. Q and Eq. pi), one can see that the one-sided quantum discord is 
not symmetric; when the measurement is performed on A, the result may be dif- 
ferent. However, there are other measures of quantum and classical correlations 
available for combined quantum systems. A symmetrical method with two-sided 
measurement over both subsystems of a bipartite system is proposed to quantify 
the classical correlation represented by the maximal classical mutual informatiorJUHl 
and the similar symmetric quantification of nonclassical correlation is definedH?3. 
A thermodynamic approach is also propose d 24 ! 25 ^ where the quantum-information 
deficit-thc difference between the total information and the extractable informa- 
tion using closed local operation and classical communication (CLOCC)-is used to 
quantify the quantum correlatiorP^. Classical and quantum correlations are char- 
acterized using measurement-induced disturbance by considering the property that 
classical states can be measured without disturbance^. Modi et al. proposed using 
the relative entropy as a distance measure of correlations-in a fashion similar to 
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the relative entropy of entanglemeirP^Lwhich provides a unified view of quantum 
and classical correlations^. A similar geometric measure of quantum correlations 
was proposed based on the Hilbert-Schmidt distance measurement. A geometric 
measurement-induced nonlocality is also defined^, which is further reformulated in 
an entropic way^. There are also other kinds of nonclassical correlation measures, 
and a recent review can be found in Ref. 1161 Usually these definitions of nonclas- 
sical correlation are not equal to each other. However, some of the measures are 
equivalent for the a special kind of Bell-diagonal statdHEH 

3. The system-environment dynamics of quantum discord 




Environment 



Fig. 1. Schematic of an open quantum system with its corresponding correlation dynamics (rep- 
resented by the wavy lines). 

As shown in fig. [I] the open quantum system S is coupled to its environment 
E, with the correlations evolving between S and E (S^E-correlation) and within 
their own systems (^-correlation and £ , -correlation) . The total Hamiltonian of the 
combined system can be written as^ 

H = H s ®lE + ls®H E + H inU (9) 

where H$ is the self-Hamiltonian of the open system S, He is the free Hamiltonian 
of the environment E, Hi nt is the interaction Hamiltonian between the system and 
the environment, and 1 represents the identity operator. 

The evolution of the closed total system (S + E), is described by the unitary 
operation U(t), where the final state becomes 

P SE {t) = U(t)p SE (0)U(t)l (10) 

The reduced density matrix of the open quantum system S becomes p$ — T^e(pse), 
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and the evolution of the quantum system can be expressed as a dynamic map 

ps(0) -► Ps (t) = Hps(0)) = Tr E {U(t)p SE U(t^}. (11) 

Fragile correlations in open quantum systems are easily destroyed by unavoidable 
noises, so knowledge of the dynamic behavior of correlations helps in designing suit- 
able protocols to protect useful resources during processing. Moreover, correlations 
can have distinctive dynamic behaviors, which can be of both practical and fun- 
damental importance. One of the most distinctive dynamic properties of quantum 
entanglement is that it may suffer from sudden death, i.e., entanglement completely 
disappears at a finite time in the system's evolutiorP. A quantum discord investiga- 
tion of open quantum systems that naturally compares other dynamic behaviors to 
entanglement obtains several interesting results. First, there is not sudden death of 
quantum discord during its dynamic phasdHEM] and almost all states have non- 
classical correlations®!. Secondly, during the dynamics of quantum discord and its 
classical correlation counterpart, there is a sudden change of behavior in their decay 
rates^. Both quantum and classical correlations are unaffected by the decoherence 
noise and the open system dynamics exhibits a sudden transition from the classical 
to the quantum decoherence regime^. Thus, in this paper, we focus on the system- 
environment dynamics of quantum discord and the fundamental influence of initial 
system-environment states on the dynamics of open quantum systems. 

As shown in fig. [I] the interaction Hamiltonian between the system and environ- 
ment leads to certain system-environment correlations that affect the states of the 
open quantum system. Maziero et al. studied the transference of correlation between 
the system and its reservoirs by considering a non-interacting two-qubit composite 
system (AB) under the influence of two independent environments (Ea and Esr^- 
By employing the Kraus operator-sum representation, the group analytically and 
numerically analyzed the dynamics of correlations of different partitions in different 
kinds of noise channels: amplitude damping, phasing damping, bit flip, bit-phase 
flip and phase flip. For the amplitude-damping channel-the dissipative interaction 
between the system and the environment that includes an exchange of energy-the 
action of the channel over one qubit is represented as 

|0>s® |0)j5 -> |0) s ® \0) E , 

\l)s <8 \0)e -> yfi^PMs <» \0)e + Vp\0)s ® \1)e, (12) 

where \0)g (|l)s) is the ground (excited) qubit state and \0)e (|1)e) represents the 
state of the environment with zero (one) excitation distributed over all its modes. 
The initial state of the whole system can be written as 

1 3 

Pabe a e b = j(l AB + )®|00) £aBb , (13) 

i=l 

where 1 00) EaEb is the vacuum state of the environments and <Xj (i — 1,2,3) are the 
three Pauli matrixes. The coefficients Cj are real numbers and the state of the system 
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Pab represents Bell diagonal states. The corresponding reduced-density operator is 
obtained by taking the partial trace of pabe a e b - 

Maziero et al. considered the correlation dynamics of AB, AEa, AEb and 
EaEJ^- Note that, due to the symmetry of the system, the density matrix of 
the partition BEb and BE a are identical to that of AEa and AEb, respectively. 
For a Werner initial state where c\ — C2 = C3 = —a (0 < a < 1), both the classical 
and quantum correlations of AB vanished in the asymptotic limit. But the entan- 
glement suffers from sudden death at a given time evolution. This phenomenon of 
entanglement sudden death between the qubits, and entanglement sudden birth be- 
tween the reservoirs, was also observed in Ref. [40j The vanishing of classical and 
quantum correlations between the parts of system AB was accompanied by the 
simultaneous creation of corresponding correlations between the reservoirs EaEb- 
Because the qubits were under the influence of the local environment, the quantum 
correlations were created between each qubit and its own reservoir, i.e. betwen AEa 
and BE B - 

The authors also investigated^ the phase-damping channel case, which describes 
the loss of quantum coherence without loss of energy. In this case, the interaction 
map for one qubit case becomes 

|0>s® \0) E -> |0) s ® |0) E , 

|l)s ® |0) B -> V^PWs ® |0) B + Vp|1) s ® |1> B . (14) 

For an initial separable system state with an independent environment, there is no 
entanglement created between the system and its own environment; i.e., AEa and 
BEb are still separable. However, decoherence does occur for the asymptotic decay 
of the off-diagonal elements of the initial state. As a result, information transference 
between the system and the environment relies on both classical and quantum 
correlations. For the phase-damping channel case, quantum correlations disappear 
in all partitions, a result different from the amplitude-damping case, where the 
quantum correlations are completely transferred from system AB to reservoir EaEb 
at p = 1. The effect of the bit-flip, bit-phase-flip, and phase-flip channels is similar to 
that of the phase-damping channel, which destroys the phase information without 
an exchange of energy. Therefore, the dynamic behavior of the correlations in the 
these three channels is the same as in the phase-damping channeP^. 

Ge et al. analyzed^ the case of two spins interacting independently with their 
own boson reservoirs with the corresponding Hamiltonian given by 

2 N 2 N 

h = £(f < + $>4,Afc) + EE^( ff ri + 4hk), (i5) 

i=l k=l i=l k=l 

Here, N — » 00. 6^.^ (i — 1, 2) is an annihilation operator of the k-th model in the i-th 
reservoir with the corresponding frequency u>k- erf, ct 2 + , (oj~) and uj represent the 
Pauli operator, raising (lowering) operator, and Zeeman splitting of the i-th spin, 
and Qij. denotes the coupling strength between the fc-th model in the i-th reservoir 
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and the corresponding spin. The group first considered an initial state with two 
excitations in the spin system 

|*o) = H0) S1 |0) S2 +/3|l) sl |l) S2 )|0) ri |0) r2) (16) 

with | ck | 2 + |/?| 2 = 1. The collective state |0) n = E[fc=i l^fc)^ indicated there was no 
excitation in reservoir n. By introducing the collective state |l) r< representing the 
case of one excitation in the reservoir, they obtained the reduced density matrix 
of p Sl s 2 , Psirn Ps 2 r 2 an d Pnr 2 - The authors then looked at two kinds of spectral 
distribution. For a flat spectral density, both quantum and classical correlations 
initially stored in the spin system decayed monotonously, only gradually transferring 
to the reservoirs. On the other hand, the Lorentz form of spectral density resulted 
in a strong non-Markovion effect. The group also found correlations between all 
bipartitions of the system and reservoir oscillations, and that both the quantum 
and classical correlations stored in the spin system also transferred to the reservoirs 
over time. During the evolution of this kind of initial input state, quantum and 
classical correlations are equal to one another in the two-spin system and the two 
reservoirs. 

In this work, another state with only one excitation in the spin system is further 
considered^, i.e. 

|*o) = («|0) S1 |1> S2 +/3|l) sl |0} S2 )|0) ri |0) r2 . (17) 

For this case, quantum correlation is no longer equal to classical correlation in the 
spin system or the reservoirs. Both the quantum and classical correlations initially 
stored in the spin system transfer to the reservoirs at the end of the system's 
evolution. 

The transference of quantum correlation between systems and reservoirs has 
also been considered in a multi-qubit systems. Man et alW^ studied a simple system 
of three entangled qubits (A, B, and C) which were coupled to three independent 
reservoirs (a, b, and c), and investigated the time evolution of discord for any two 
qubits. The Hamiltonian for this total system can be expressed as the sum of three 
independent spin-boson models 

H = H Aa + H Bb + Hqc (18) 

where H Xx = u q o\o x +J2 3 -^] x ) x j + J2j(9j<^x x j + 9] (J x x t)' Xx e {Aa,Bb, Cc}, 
which is similar to the Hamiltonian of Eq. ( [15] ). The qubits A, B, and C are initially 
prepared in the extended W-like state and read as 

PABC ABc{W(i))\, (19) 

where \W(0)) ABC = a|001) +/3|010) + 7 |100) with |a| 2 + |/3| 2 + |7| 2 = 1, and 1 ABC is 
the identity matrix. The three reservoirs a, 6, and c are initially in the vacuum state 
1 000) a b c , and the initial total system is in the product state PABCabdfy = PABcify® 
|000) o (, c (000|. The reduced matrix operator of any two qubits, e.g. qubits AB, is 
obtained by tracing the total evolved state PABCabc(t) over the reservoirs abc and 
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the qubit C; the other reduced density matrixes can be obtained in a similar fashion. 
Man et oiP^ first considered a Markovian evolution case with a — f3 — 7 = 1/ \/3- 
When p = 1 (W-state), the group found that both the entanglement and discord 
of qubits AB decayed monotonously to zero and transferred to the corresponding 
reservoirs ab. On the other hand, for p = 0.8 (mixed state) and p = 0.5 (separated 
state), the quantum discord initially in the qubits AB transferred to the reservoirs 
ab but with a sudden change in the behaviors of their dynamics rates. The group 
showed that while entanglement can suffer from sudden death and sudden birth, 
quantum discord always evolves asymptotically. The authors further considered^ 
a model with two remote nodes, each of which consisted of two qubits coupled to 
a common reservoir, i.e., qubits AC in reservoir a and qubits BD in b. The initial 
state of the total system was prepared to be pACaBDb{0) = Pab(0) <8> Pcd(0) <8> 
Pab (0), where p AB = (1 - p)/Al AB + p|*)ab(*| with \V} AB = a\00) + j3\ll), 
Pcd = |00)cd(00|, and p ab — |00) a {,(00|. As a result, the two nodes were initially 
correlated by the correlations between the qubits A and B. Both Markovian and 
non-Markovian cases were investigated, with the group finding that the quantum 
discord of AB, CD, and ab all tended towards long-term steady values, while the 
discord of ab was always larger than that of AB and CD, and the entanglement 
of AB and CD eventually decayed to zero. Similar to the first model, there were 
sudden changes in behavior in the quantum discord of the qubits and the reservoirs 
when the initial states of AB were prepared in a mixed system. 

In the above discussion, the environments were modeled as the quantum system. 
However, what is the result when the environment is not quantum but classical? Lo 
Franco et al. considered^ a pair of noninteracting qubits, each locally coupled to 
a classical random external field. The external field was equal for both qubits and 
was unaffected by the coupled qubit. The amplitude of the external field was fixed, 
but the phase of each mode was classically random and equal to either or 7r with a 
probability of 1/2 (the case for which p 7^ 1/2 was also considered by the authors). 
For each external field with phase or 7r, the qubits underwent a unitary evolution. 
The dynamical map for the initial two qubits pab can be written as (whenp = 1/2) 

PAB® = 4 E U A (t)Uf P AB{Q)U^Uf\ (20) 

with Uj(t) — e r lE -i t l h _ Here, (S = A,B) represents the time evolution operator, 
where Hj = ihg(a + e^ l ^ j — cr_e 1 ^). The time evolution operator Uj(t) can be 
expressed in the basis {|0), |1)} as 

T tSm-{ cos ^ e-*+* em.(gt)\ . . 

U i [ ) ~ \ -e^ sm(gt) cos(gt) J ' [ ' 

where j = 1, 2 with <f>i = and 4>2 = tt- 

When the initial state was prepared in the Bell diagonal state of pab = 
A+|^ + )<^ + |+A_|^_)<V-|,with A+ =0.9 and A_ =0.1 (|V'±) = l/%/2(|01) ± |10))), 
the authors observed the collapse and revival of entanglement, quantum discord, 
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and classical correlation, with all the correlations measured in the unified relative 
entropy forrrPsl. The group determined that classical correlation oscillates when 
quantum discord remains constant, and vice versa. Thus, classical correlation and 
quantum discord are equal to each other at certain times in the system's evolution, 
which is the same as the case of a sudden transition between classical and quantum 
decoherence^. 

It is surprising that the classical environment has no back actions and does not 
store the quantum information of the two qubits. Tasking into account the corre- 
lations in the quantum-classical state of the system-reservoir system, the authors 
suggested that the recorded operations applied to the qubits in the classical envi- 
ronment could play an important role in retrieving the quantum correlations^. 
The group also discussed a class of global system-environment evolutions pos- 
sessing a lack of back-action. For this class, any reduced system dynamics ob- 
tained with a quantum environment can also be obtained by modeling the classical 
environment!^. In order to interpret the revival of correlations, the authors general- 
ized the non-Markovian quantifier introduced in Ref. 33] by using a time-dependent 
non-Markovianity quantifier 

lRn(t) = ^ dRn[ P Mt')] w _ Ai?n(;); (22) 

J to dt 

Here, Rn represents the relative entropy of entanglement)^, with ARn(t) — 
Rn(t) — Rn(to) where to < t. The authors observed that I Rn (t) monotonically 
increased in the same time regions when revivals of entanglement and quantum 
discord occurred^. 

The initial correlations between the system and environment can fundamentally 
influence the dynamics of the system. One of the most famous examples of this is 
the completely positive map for describing the dynamics of open quantum systems^. 
If the dynamical map Eq. (Illl) can be expressed as 



ps(t) = *(ps(0)) = J2 a i C iPs{0)Cj, ai >0Vj (23) 

i 

where Ci is the corresponding matrices and J^i CiC\ = 1, it is a completely positive 
map^, and thus can described the dynamics of the system when the initial state 
between the system and its environment is a product stated In a recent work, 
Rodriguez et al. extended^ this result and demonstrated that, if the initial state 
between the system and environment is a classically correlated state with zero dis- 
cord, the dynamics of the open quantum system are a completely positive map. 
The system-environment state with vanishing quantum discord can be represented 
as pse = ^2iPiMf ® PeHi where pe\i are density matrices for the environment and 
{Mf} is a complete set of orthogonal projectors on the system S with Pi > and 
^2iPi = 1. For such initial system-environment states, Rodriguez et al. proved^ 
that, under the unitary evolution of the total systems, the final state of the system 



can be expressed in the form of Eq. ( 23 ) , which implies a completely dynamical 
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map. By introducing a kind of special-linear state with the property of being of 
full measure in the set of mixed bipartite states, Shabani and Lidar demonstrated^ 



that a quantum dynamical process with the form of Eq. ( 11 ) is always a linear Her 



mitian map for arbitrary initial system-environment states (pse(0)). They further 
proved that the vanishing quantum discord of pse(0) is a necessary and sufficient 
condition to induce a completely positive map. 

Coupling between the system and the environment leads to decoherence and 
changes the entropy of the system, and von Neumann entropy S is typically used to 
quantify this decoherence. Thus, the rate of decoherence is related to the entropy 
rate of the system, which is defined as ^ . Recently, a general result for determin- 
ing a system's entropy rate has been demonstrated^. The sufficient and necessary 
condition for a zero system entropy rate at a time r under any system-reservoir 
coupling is expressed as 



dt 



V H & [p s 



l E ,p SE ] 



0. 



A system-environment state with the property of [p s 
lazy state, which has the form of p SE = ^^Mf ® 1 



g) l E ,p SE ] 
E p SE Mf ( 



(24) 

= is called a 
1 E . Note that 



p s = YljPjMj and {Mj } are orthonormal projectors. As a result, the lazy state 
is a classical state with vanishing quantum discord^. The authors further provided 
a universal bound on the rate of decoherence for a system-environment interaction 
Hint of arbitrary strength, which is expressed as 



dt 



< \\Hi, 



[\n(p s ) ® 1 E , p 



SE1 



(25) 



Due to the fact that almost all quantum states have nonclassical correlations^, 
one may conclude that the rate of system entropy will be always high. However, 
in another recent worlP3 the authors demonstrated that this is not the case. The 
authors showed that the probability for obtaining a high entropy rate for a randomly 
chosen state pse is very small, i.e. 

d 



Pr 



PSE 



dt 



S > Hi 1 



<S, 



(26) 



where e = 2- 1 /2(iog 2 d E -3io g2 d s -4) and s = 2e -d s /i6 with ds and d]s representing 
the dimensions of the system (S) and the environment (E), respectively. For a 
sufficiently large environment (log 2 ds > 31og 2 ds) and not too small of a system 
(log 2 ds > 2), the authors found that the entropy rate of S in a vast majority 
of bipartite states can only decay at a vanishing rate. For small systems where 
log 2 ds < 2, a similar bound exists when log 2 dE > (9/2) log 2 dj^l. 



4. External control effects on open quantum systems 

One can greatly influence the dynamics of correlations by implementing opera- 
tional control on an open quantum systems. In the work of Francica et al, the 
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group considered^ quantum Zeno and anti-Zeno effects on quantum and classical 
correlations of two noninteracting qubits coupled to a commonly structured reser- 
voir, e.g. two atoms coupled to a common zero-temperature bosonic reservoir. The 
authors performed nonselective measurements on the system AB at time interval 
T to observe the quantum Zeno effect on entanglement, and found two characteris- 
tics: 1) one of the possible measurement outcomes is a projection onto the collective 
ground state \ipo)s = |0)a|0)b, and 2) the measurement cannot distinguish between 
the excited-states \ipi) s = |1)a|0)b and \ip 2 ) s = lO^ll)^- 

The measurements can thus be described by the two following projectors 

Mo = \ip )s(M ® 1b, 

Mi = (|Vi>s(Vi| + \1*)s{\H) ® 1b, (27) 

In this way, the group separate results for the different cases of off-resonant interac- 
tion between the qubits and the cavity mode and the resonant casein The authors 
observed that, under certain conditions, the same measurements that greatly sup- 
pressed entanglement loss could instead lead to a much faster disentanglement, 
showing the anti-Zeno effectlHEH the group further found that the dynamics of 
quantum discord and classical correlation are qualitatively similar to the dynam- 
ics of entanglement, and a series of oscillations between the Zeno and anti-Zeno 
effects can occur as a function of the time delay between successive measurements. 
Studying the fact that classical correlation can be protected by a sufficiently high 
measurement frequency that can also display the Zeno and anti-Zeno oscillation 
could lead to a deeper understanding of classicality. 

Bang-bang control technology, as a dynamical decoupling strategy, has been used 
to prevent quantum decoherence^. Xu et al. applied the bang-bang control pulses 
to two noninteracting two-level atoms in an independent environment to investigate 
its influences on the dynamics of quantum correlations-^!. When adding the bang- 
bang pulses, the Hamiltonian of one atom interacting with its own reservoir can be 
expressed as 

H = H + H in t + H p , (28) 

with Hq — ^j-<Tz + coa^a and Hi nt = g(a + a + er_a^), where a and a) denote the 
annihilation and creation operators for the cavity held. The Hamiltonian for a train 
of identical pulses of duration r is H p , expressed as^ 

oo 

H p = Va z Y^Oit-T- n(T + r))6»((n + 1)(T + r) - t), (29) 

n=0 

where V is the amplitude of the control field and T is the time interval between 
two consecutive pulses. In the authors' case, V was set to 7r/2r, which implies the 
7r-pulse only. The two atoms were initially prepared to be the Werner-like states 

PAB = xm^\ + 1 ^i, (30) 
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where |$) = o;|0)a|1) b +/3|1)a|0)b with \a\ 2 + |/3| 2 = 1 and x represents the purity 
of the initial states. The cavity field was prepared in the Fock state or the thermal 
state. The group determined that both entanglement and quantum discord can 
be enhanced by the application of bang-bang control, and the increased amount 
is larger with shorter time intervals (T) between the control pulses 56 . When the 
cavity field was initially in the vacuum state, at time t2N = 2N(T + r), quantum 
discord was found to recover to the initial value for the detuning 5 — (S = ujq — uj) . 
When 5^0, the quantum discord fluctuated at time t2N with period 2ir/5 and the 
oscillate amplitude was smaller with the shorter interval. The maximal revivals of 
quantum discord in the thermal state case decreased slightly with time, and the 
decrease became even slower with shorter control pulses. 

Li et al. investigated^ the effects of practical feedback control on the dynamics 
of entanglement and geometric quantum discorcP^ of two identical atoms (A and B) 
resonantly coupled to a single-mode cavity by considering the Markovian feedback^. 
In this system, the control Hamiltonian is H c = I(t)F, where I(t) is the signal 
outputted by the homodyne detection and F is the feedback Hamiltonian. Taking 
the efficiency of the detection, denoted by r), into consideration, the modified master 
equation for the state of the two atoms becomes^! 

f = -i[H + \{&F + FG), p] + D[G - iF]p + D[^^F]p. (31) 

where H = Q(<7^ + <J^) represents the driving of the laser and D[X]p — 
XpX' — (X'Xp + pX'X)/2 represents evolution due to coupling with the oper- 
ator X between the system and the environment. Note that G = r(cr^ + cr?), 
which represents the jump operator (r was set to be 1 in this work). If both atoms 
are initially in the excited state (|11)ab), entanglement first increases and then 
decreases with feedback control. The phenomenon of entanglement sudden death 
occurs when 77 < 1, and the duration of the non- entanglement period increases as 77 
decreases. Geometric quantum discord displays a similar tend to entanglement, but 
without sudden death. When both atoms are initially in the ground state (|00)ab), 
the noise of the detector (77 < 1) can trigger the detector and create entanglement 
and quantum discord between the two atoms. For the initial entangled input state 
(111)^3 — |00)^s)/v2, the authors found that there was also a sudden change of 
entanglement with the perfect detection due to the discontinuous maximum value 
in calculating concurrence®!. 



5. Experimental investigation of the dynamics of quantum discord 
in two-qubit systems 

Several distinctive properties of the dynamics of quantum and classical correlations 
under decoherence were experimentally investigated by Xu et aZP^ using a phase- 
damping channel simulated by birefringent quartz plates as the noisy environment. 
The interaction quantum map was similar to that of Eq. ( 14 ). The information carri- 
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ers were encoded as photon polarizations with the coupling between the photon's po- 
larization and frequency modes occurring in the birefringent quartz plates. This cou- 
pling led to a phase-damping effect through tracing the frequency freedoms. By con- 
structing unbalanced Mach-Zehnder devices, the group prepared different kinds of 
Bell-diagonal states, i.e. p A B = ci|$ + )($ + |+c2|$~)($~|+c 3 |* + )(* + |+c4|<I'-)(1'-| 
with |$±) = 1/V2(|00> ± |11)) and = l/v^flOl) ± |10)). The minimum con- 
ditional entropy of Eq. ^ was obtained by scanning the measurement on photon 
B. Quantum discord and classical correlation were obtained from the final states 
reconstructed by quantum state tomography^. For an initial state with c 2 = 3/4 
and C4 = 1/4 (ci=0 and C3=0), there is a sudden transition from the classical 
decoherence to the quantum decoherence regime^, as shown in figure 2 (the tran- 
sition point is represented by the crossover of the four panes with different colors). 
The figure also shows the non-entanglement quantum correlation, which exhibits 
a sudden change in the decay rate^-l Quantum discord monotonically decays after 
the sudden transition point, while entanglement suffers from sudden death after 
a certain period. The experimental results are in good agreement with theoretical 
predictions. For a different initial state with c\ = 0.09, c 2 = 0.09, c 3 = 0.81 and 
c 4 = 0.01, quantum correlation is larger than classical correlation as the system 
evolves, disproving an early conjecture that classical correlation is always larger 
than quantum correlation^. 




J i I i I i I i I i I i L 



50 100 150 200 250 300 (\) 
Time Evolution 

Fig. 2. The correlation dynamics of the input state pab = § I* - }{®~ I + 2 1* - } (* - 1- The green 
upward-pointing triangles, black squares, red dots, blue stars, and magenta downward-pointing 
triangles represent the experimental results of T (mutual quantum information), C (classical cor- 
relation), Q (quantum correlation), En (entanglement of formation) and Rn (relative entropy of 
entanglement), respectively, with the green solid line, black solid line, red dashed line, blue dotted 
line and magenta dotted line representing the corresponding respective theoretical predictions. 
Non-entanglement quantum correlation (T>) is compared to C in the inset; purple dots represent 
the experimental results of T> and the purple dotted line is the corresponding theoretical prediction. 
The x-axis represents the total thickness of the quartz plates, with Ao = 0.78 (im corresponding 
to the time evolution of the photons. This figure is reproduced from Xu et al. Nat. Comm. 1, 7 
(2010). 
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The non-Markovian dynamics of classical and quantum correlations was fur- 
ther investigated in an all-optical setup by Xu et alP^. In their experiment, the 
group simulated a non-Markovian channel using a Fabry-Perot cavity followed by 
quartz plates acting on one of the two photons. The Fabry-Perot cavity filtered the 
frequency distribution of the photon, leading to a refocusing effect of the relative 
phase in the dephasing environment and inducing the non-Markovian effect. The 
sudden transition from the classical to the quantum decoherence regime occurred 
at the beginning of the evolution, similar to the case in fig. 2. Under such a non- 
Markovian environment, the quantum correlation (measured by the relative entropy 
of discord^ and calculated from the reconstructed density matrixes) revived from 
near zero before decaying again in the subsequent evolution. However, there was 
not revival of the classical correlation, which remained constant over the same pe- 
riod. This phenomenon is explained by the weak non-Markovian effect. The authors 
theoretically showed that both classical and quantum correlations are revived with 
a narrower frequency width (i.e. stronger non-Markovian effect). The group further 
implemented an a x operation on the photon under decoherence and investigated the 
corresponding correlation dynamics. This process found a sudden transition from 
the quantum to the classical revival regime, as well as correlation echoes. 

Experimental investigations of quantum discord and classical correlation under 
decoherence have also been implemented using NMR (nuclear magnetic resonance) 
systems. Because of the high-temperature expansion, a typical NMR-system den- 
sity matrix for the effect two-qubit case can be written as p = |l + eAp, where 
1 is the identity matrix, e is the thermal energy, and Ap represents the deviation 
density matrbP^ that can be reconstructed by quantum state tomography. Under 
the action of a global environment, Soares-Pinto et alW* investigated the dynamics 
of correlations in a NMR quadrupolar system. The symmetric classical and quan- 
tum correlations' 7 ! 21 ! , which were computed from an experimentally reconstructed 
deviation-density matrix, were found to both decay exponentially. These experimen- 
tal results agree well with the theoretical predictions. The authors further showed 
that, depending on the initial states, the classical correlation could be larger or 
smaller than the quantum correlation during the system evolution. 

The sudden change in behavior in correlation dynamics has been further demon- 
strated in the NMR systems at room temperature^. After preparing different kinds 
of initial Bell-diagonal states with specific relations between their components, the 
authors left the systems to evolve freely under decoherence. The environment of 
such an NMR system can be modeled by two independent quantum channels: the 
phase-damping and amplitude-damping channels. The authors observed a sudden 
change phenomenon for quantum discord and its classical counterpart due to the 
presence of phase-damping noise. On the other hand, the amplitude-damping noise 
led to a small decay of the classical correlation compared to a constant classical cor- 
relation after the sudden change point. Note also that the two nuclear spins finally 
relaxed to the Gibbs states. 

Experimental investigation of quantum and classical correlations usually requires 
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Fig. 3. Schematic setup for witnessing nonclassical correlations in a two-qubit system AB under 
dccoherence. This figure is reproduced from Auccaise et al. Phys. Rev. Lett. 107, 070501 (2011). 



full quantum state tomography and their resulting numerical calculations. However, 
the concept of witnessing the quantumness of correlations recently received a great 
deal of attention. For a two-qubit states with the form of pab = (1 AB +Xh=i {a-i^t® 



-bil A ®<T?- 



■Ci<jf (E><Ji )) / '4, Maziero and SerraP^ obtained a sufficient condition 



for the states to be only classically correlated 



W 



t f \<Oi 

i—1 j=i-j-X 



,(Oj) P \ = o, 



(32) 



UA^B[ReM)P^i))]UA-,B. In this 
A(B), 



In this equation, Oi = af ® af , with i = 1,2, 3, and O4 
with z,w £ 3? 3 and a £ {01, 02, 03}. Note that Eq. (32 1 is also shown to be a 
necessary condition for the Bell-diagonal states. The value of (erf 
can be changed to (af ® l B ) 7i with 7,, 

equation, R ei (0i) = R A (9i) ® Rf.(6i) with R2^ a> [di) representing a local rotation 
with angle Q{ around direction e; on qubit A (£?), 6*i = 0, 6*2 = 6*3 = 7r/2, e2 = y, 
e3 = z, and Ua^b represents the controlled-NOT gate between qubits A and B with 
A as the control qubit. Auccaise et al. experimentally demonstrated^! such a witness 
in a NMR system at room temperature by preparing different kinds of Bell-diagonal 
states (in this case (O4) = 0) and performing a corresponding witness circuit, and 
the schematic setup is shown in fig. [3] (where the evolution process is added in 
the case considering the witness dynamics). For each initial state, the experimental 
process was carried out three times to obtain the values (ct] 4 ^; for determining the 



witness in Eq. ( 32 ) . Quantum state tomography and numerical optimization were 
used to calculate correlation values for comparison, and they agreed well with the 
witness results. The authors further investigated the decoherence dynamics of the 
witness with a freely evolving initial state. The witness tended to be 0, with the 
two-spin state evolving to a state with only classical correlation. 

Experimental investigation of quantum discord has now been extended to solid 
systems. In a recent experiment, a system of electron and nuclear spin in a phos- 
phorous donor in silicon was observed by generating a series of separable thermal 
states with different anisotropic parameters^!. The authors observed nonzero quan- 
tum discord and found that there were sudden change behaviors of the quantum 
discord during tuning of the anisotropic parameters. There are also other exper- 
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imental works concerning total quantum correlations; e.g. a modified correlation 
witnes£3Z3 and the space of quantum correlations^. These types of systems could 
provide an intriguing perspective on the dynamics of quantum discord, and help to 
better determine different experimental technologies. 

6. Conclusion 

Past investigations of quantum discord in open quantum systems have helped to 
disclose the particular dynamical behaviors of discord, as well as applications of 
both fundamental and practical importance. Nevertheless, observations of Marko- 
vian zero-discord classicalitjP^, and the unusual response of quantum discord to 
temperature^, still show that further investigations on the quantifications of quan- 
tum and classical correlations are required. Another arresting aspect is that quan- 
tum discord can be created under local operations' 7 ' 6 ! 77 ! 7 ' 8 !, which is different from 
that of entanglement^. Gessner et al. suggested^ that quantum discord should be 
interpreted as the local measure of quantumness, not necessarily as the number of 
quantum correlations. As a result, the problem of classifying quantum and classical 
correlations is still open despite the existence of the many different ways to measure 
nonclassical correlation introduced in Part 2. 

From the above theoretical and experimental investigations, we can see that the 
time evolutions of both quantum discord and classical correlation in open quantum 
systems can be deduced from the time evolution of the initial states by using quan- 
tum state tomography and optimization calculations (or the corresponding witness 
circuits) . Inspired by the factorization law of entanglement evolution in noisy quan- 
tum channels' 80 ! 81 ! 82 !, we expect that a similar simple relationship exists to char- 
acterize quantum correlation in a noisy environment. It has already been shown 
that any quantum correlation should be related to the amount of entanglement!^. 
Therefore, suitable quantification of nonclassical correlations is needed to meet this 
mission. 

More experimental works must be done to further understand quantum discord, 
especially to clarify the quantum advantages in quantum information processing 
and to demonstrate their fundamental applications. However, increasing interest in 
this type of work will certainly lead to more distinctive discoveries in this field. 
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